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It is proposed to compute matrix elements for the (unphysical) K^n^ transition to deter- 

mine the next-to-leading order low energy constants of the weak chiral Lagrangian. This allows 
us to evaluate K — * {nn)j=Q decay amplitudes at this level of precision. This approach has sev- 
eral significant advantages over the use of A' nn transitions, most notably the elimination of 
s-channel disconnected diagrams and the use of fewer inversions. 
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1. Introduction 

Nonlep tonic kaon decays present a number of major cliallenges. Direct CP violation in kaon 
decays has been confirmed by the non-zero measurements of e'/e [1,2], but a quantitative theoret- 
ical understanding, even within the standard model, is still lacking. The AI =1/2 rule also remains 
a longstanding puzzle. Whilst these are weak decays, the relevant energy scale is of the order of a 
few hundred MeV at which non-perturbative QCD effects are significant making the evaluation of 
the amplitudes difficult. 

The lattice formulation of QCD together with large scale numerical simulations provide the 
opportunity to compute the non-perturbative effects from first principles. However, with the com- 
puting resources which are currently available it is not possible to perform the simulations directly 
at the physical values of m„ and rud- Although the situation is rapidly improving, it is still unlikely 
that computations with pion masses of about 140 MeV will be performed in the near future. It is 
proving particularly useful in general, and for kaon decays in particular, to use chiral perturbation 
theory (^PT) to obtain physical results from simulations with heavier values of mu,d [3]. In this 
approach, the QCD operators in the effective weak Hamiltonian are written in terms of a sum of 
operators composed of meson fields in the same chiral representation. A priori, the coefficients of 
the operators in the resulting Weak Chiral Lagrangian, i.e. the low energy constants (LECs), are 
unknown. The aim is to determine the LECs from the mass and momentum dependence of matrix 
elements computed in lattice simulations. This is well understood for A/ = 3/2 K ^ TlTl decays in 
which the two-pion state has isospin 7 = 2. In this paper we focus on A/ = 1/2 decays in which the 
two pions have 7 = 

At leading order in ;tPT, there are only a small number of LECs. For example, for the A5 = 1 
operator in the (8,1) representation there are 2 LECs, ttj ' ' and , in terms of which the 
K 7171, K ^ 71 and K vacuum matrix elements are: 

(;r+;r-|^(«'i)|^) = ^(ml-m^Jaf^ (1.1) 

(;r+|^(«'i)|i^«) = ^(a^)-a(^'i)), (1.2) 

(0|^(«'i)|^») = ^(m|-m^,)a(^'^ (1.3) 

Thus by computing the K 7i and K vacuum matrix elements, a\ 2 can be determined and 
the K ^ 7171 matrix elements can be evaluated at leading order. Such calculations were performed 
in 2001 with quenched ensembles and with meson masses above about 600 MeV [4,5]. While 
these calculations demonstrated the feasibihty of the procedure, it is clear that, in addition to using 
dynamical quarks it is necessary to perform simulations at fighter masses and to go beyond the 
leading order in the chiral expansion in order to understand the A/ = 1/2 rule and the value of 
e'/e. The increase in computer power since 2001 and Next-to-Leading Order (NLO) calculations 
in ;^PT [6-10], make it possible to contemplate such an endeavour. At NLO in ;tPT, the number 
of LECs to be determined grows and it is not possible to determine them all from K ^ 71 and 
K vacuum matrix elements alone. 

One possible approach to the determination of the NLO LECs is to compute K ^ 7i7i matrix 
elements. However, as explained in the following section, there are difficulties in studying two- 
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pion states with isospin 1=0. In this talk, we present an alternative procedure for determining all 
the necessary LECs at NLO, based on the evaluation of K^K'' transitions (see sec.||). A 

detailed analysis of this proposal will be presented in a paper which is currently in preparation. 

2. Difficulties in the Evaluation of ^ — > {kk)i=o Matrix Elements 

In this section we discuss the evaluation of K ^ nn matrix elements and the difficulty en- 
countered by the presence of ^-channel disconnected contractions. For illustration, consider the 
process K*^ n^n^, studied in ;^PT in refs. [7, 8]. Using the operator product expansion, the 
amplitude can be written in terms of matrix elements of the A5 = 1 weak effective Hamiltonian 
which contains 10 four-quark operators 2, (/ = 1-10): 

G 

{7t7t\JifAs=i\K) = ^£ V^^^c^At) {7im\K)^ , (2.1) 

where V^j^j^^ are appropriate combinations of CKM matrix elements, c, (/x) are Wilson coefficients 
and /I is a renormalization scale. To demonstrate the ideas let us consider one of the operators: 

Qi = SaY^i^ - 75)da Mfo7^(l " /s)"/.- (2-2) 
The correlation function from which the K ^ nn transition is determined is 

CK^7in{tK,to,t7i) = (0 \^27t{h) Qiito)ffK{tK)\0) (2.3) 

« {0\^27i\7i7i){7i7i\Qi\K){K\ffK\0)e-^'"'^'''-'°'>e-'^'''^'''-'''\ (2.4) 

where in the last line we assume that the time intervals tjc — to and to — tx are sufficiently large and 
positive so that contributions from heavier states can be neglected. &k and i^27i are interpolating 
operators for the kaon and two-pion states respectively. By computing Ck^uji one can check the 



validity of expected asymptotic dependence on the time intervals in eq. ( |2.4j ) and extract the matrix 
element {nn\Q\\K). 

The quark flow diagrams which contribute to the correlation function are sketched in Fig.|l|. 
The number of quark propagators, and hence inversions of the Dirac operator, required to evaluate 
these diagrams depends on the strategy as to which of the times tx, to and tj^ are chosen for the 
sources. If we choose to allow to vary, then we have to solve the problem of inserting the 
quark propagator connecting 71+ and tt" in diagram (A,B,C) in Fig.[T[ This requires as many quark 
propagators as the number of time slices f^r. If instead we allow to to vary, the insertion of the loops 
in diagrams (A) and (C) require a large number of inversions. If we fix both to and tji then we are 
not able to check the time behaviour of the two-pion state. Even if it were possible to perform the 
large number of inversions which are required to evaluate the diagrams, from experience we expect 
that a very large number of configurations would be required to evaluate the disconnected diagram 
(C). Indeed one might try evaluating the diagrams using stochastic all-to-all propagators, but this 
may be noisier. 

We stress that the difficulties described above are technical rather than fundamental. They are 
nevertheless delaying the evaluation of matrix elements with two-pion states in the 1=0 channel (in 
the 1=2 channel only diagram (D) needs to be evaluated, which is relatively easy). 
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Figure 1: Quark contraction diagrams for K nn matrix elements. The grey circle represents the insertion 
of a four-quark operator. 
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Figure 2: Quark contraction diagrams for the propagator of the n^n state. 



In order to extract the matrix element {nn\Q\ \K) from Ck^tik we need to divide by the overlap 
factors, including {0\tffjc7i\^^) (see eq. (2.4)). The evaluation of {0\i&jc7i\^^) requires the computa- 



tion of the propagator of the 1=0 two-pion state, i.e. the evaluation of the diagrams in Fig.|^ This is 
at least as challenging a task as that of the Ck^hh, with similar issues to the ones described above. 



?>. K^'n ^ n Matrix Elements 



We start by stressing that our aim in this paper is limited to the determination of the LECs at 
NLO in the chiral expansion, from which we evaluate the K ^ nn decay amplitudes at the same 
level of precision. The key point is that it is possible to compute the LECs using any appropriate 
external states; in particular it is not necessary to use two-pion states with 7 = 0. Our proposal 
is to determine the LECs for A5 = 1 A/ = 1/2 operators by computing matrix elements for the 
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Figure 3: Quark contraction diagrams for Kn n transitions. The grey circle represents the insertion of a 
four-quark operator. 

(unphysical) transition K^^n^ —>■ as well as ic' ^ 71 and K vacuum matrix elements. In this 
section we explain the advantages of this proposal; a detailed demonstration that all the necessary 
LECs can be determined from these transitions for a reasonable set of kinematic parameters will be 
presented in a forthcoming publication, ^"'tt^ is an / = 3/2 highest weight state, and, as we shall 
see below, this leads to a number of important simplifications (analogous to those present when 
studying the I = 2nn state). 

The correlation function corresponding to Kn n transitions is 

CKOn-^K-itK,to,tn) = {0\ffnAt^)Ql{to)^n-{tK)^K{tK)\0) (3.1) 

where ffj^- and are interpolating operators which can create or annihilate a n~ respectively. 
The quark contraction diagrams corresponding to this process are shown in Fig.|3|; they have com- 
mon topologies with those in Fig.|I|. The difference is that the final-state 71+ in Figjljis now crossed 
into the initial state where it is a tt^. Thus, the quark line connecting 71+ to in diagrams (A,B,C) 
of Fig.[l] now connects 71^ to propagating from Ik to tj^ in Fig.|[ We can therefore vary Ik using 
a small number of quark propagators, and hence check the asymptotic exponential behavior of the 
Kn state. 

What is more important is the absence of diagram such as (C) in Fig.[l|. The i^-channel discon- 
nected diagram in Fig|I](C) now becomes the f-channel disconnected diagram in Fig.^C). Such 
^-channel disconnected contributions appear, for example, when studying K ^ nn matrix elements 
with an / = 2 7171 state. There have been several calculations of I = 2 nn states [11, 12] and this 
did not pose any problem. 

A related advantage is that one can avoid the vacuum subtraction. For an / = final state at 
rest, there is mixing with the vacuum state which must be subtracted ^ . In this case, the correlation 

'One can avoid the vacuum subtraction by considering the K7t state at non-zero momentum 
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Figure 4: Quark contraction diagrams for the propagator of the K^K state 



function in eq. (2.3) contains the vacuum contribution, 

{OWin\^){0\Qy\K''^){K^^\&K\0)e-'^'''^'o-tK) (3_3) 



which must be subtracted in order to get Eq.(2.4). There is no such subtraction necessary in the 
Kn — > 71 correlation function. 

Note also that the overlap factor (K^n^\GK^n- 10) can be calculated much more easily. Fig.^ 
shows the quark contraction diagrams. These are the same diagrams which appear in the I = 2 kk 
propagator which is calculated without any difficulty. 

When evaluating K^k'' —>■ matrix elements, we need to decompose the operators Qi into 
their A/ = 1/2 and AI = 3/2 components. InK ^ nn transitions, the isospin transfer, AI is deter- 
mined by isospin of the final state. A/ = 1/2 if Ij^n = and A/ = 3/2 if Ij^ti = 2. For the K^n^ 
transition, the isospin of the external states is already fixed, I^o^- = 3/2, = 1, both A/ = 1/2 
and AI = 3/2 transitions are allowed and therefore a basis of operators with fixed isospin should 
be used (see for example Appendix A of ref. [5]). 

Since, for both K Tin and Kn —>■ n transitions we have a two particle state, the finite- volume 
effects are not exponentially small and should be taken into account [14]. Just as for the (7171)7=2 
state, the lowest energy K^K'' state(s) can be isolated and the standard techniques [14-16] can be 
applied. Energy, of course, must be injected at the operator (which is also the case in practice in 
A' — > 7171 transitions). 

The observation that the initial ^"71^ state is one of highest weight also means that one can also 
envisage performing the Kn —>■ n calculations in the (non-unitary) partially quenched QCD, which 
is not the case for the standard extraction ofA/=l/2A'^7r7r matrix elements [13]. One might also 
try to improve the precision of the determination of the LECs by extending the kinematical reach of 
the calculations using partial twisting, which again is not possible for K —>■ (7171)7=0 transitions [17, 
18]. 



4. Concluding Remarks 

In this paper we have presented a suggestion for the determination of the LECs necessary 
for the evaluation of the K —>■ (7171)7=0 decay amplitudes at NLO in the chiral expansion. Our 
approach is based on the computation of K^n^ —>■ matrix elements (combined with K n and 
K —>■ vacuum matrix elements). We will demonstrate explicitly that all the LECs can be determined 
in this way in a forthcoming paper. We then have to develop effective strategies for implementing 
these ideas in numerical simulations and to investigate how precisely the LECs can be determined. 
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In spite of the relative simplicity of our approach, there will be technical issues in its imple- 
mentation; including the subtraction of power divergences. As is well known, some of the A5 = 1 
four-quark operators mix with lower dimensional ones, leading to ultra-violet divergences propor- 
tional to inverse powers of the lattice spacing. The numerical subtraction of the divergences at 
leading order in the chiral expansion was performed in the quenched studies of refs. [4, 5] and it 
remains to be seen how accurately this can be done in dynamical simulations at NLO in the chiral 
expansion. 
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